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ABSTRACT 

It is known that Neutron Stars may be converted into more compact Strange Stars 
(SS) on capture/formation of a "seed" of strange matter. It is also known that the 
binding energy of the nascent hot SS is likely to be radiated as v — us so that an 
electromagnetic pair fireball (FB) may be created by neutrino annihilation. But we 
show here that, a General Relativistic treatment of the problem may lead to a FB 
energy (Qfb ~ 10 53 erg) which could be higher by as much as four orders than what 
was estimated previously. Further since the baryonic mass of the thin crust of a strange 
star is negligible, this FB will generate an extremely relativistic blast wave. Thus this 
process may be one of the viable routes for the genesis of hitherto unexplained cosmic 
Gamma Ray Bursts. 
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1 INTRODUCTION 

It is now clear that a large number of Gamma Ray Bursts 
(GRBs) involve emission of 7-rays as large as Q 1 ~ 10 52 — 
10 54 erg under condition of isotropy (Kulkarni et al. 1999). 
However, many GRB afterglow light curves decline more 
rapidly than what is expected in simple isotropic fireball 
models, and it is plausible that such afterglows are actually 
beamed (Kulkarni et al. 1999, Huang et al. 2000). In such 
cases, the total (electromagnetic) energy budget could be 
substantially reduced (often by a factor ~ 100). Nonethe- 
less, many afterglows do not fade with such unusual rapid- 
ity and the decay of the light curve may be smoothly fitted 
by a single power on the time scale of months. Such GRBs 
are likely to be more or less spherical events, and thus, for 
GRB 971214, we may indeed have Q 1 » 3 x 10 53 erg (Dal 
Fiume et al 2000). There have been many attempts to ex- 
plain the energy budget and origin of powerful GRBs with 
especial emphasis on the beamed (Paczynski 1998, Popham, 
Woosley & Frier 1999, Ruffert et al. 1997). In the follow- 
ing, we put forward an alternative model to explain the ori- 
gin of GRBs with a value of Q 1 ~ 10 53 erg for which the 
problem of baryonic contamination is minimal because we 
are considering the possibility of conversion of an isolated 
or mass accreting massive neutron star (NS) into a Stange 
Star (SS). It might appear at first sight that the model is 
not new at all because starting with Alcock, Farhi & Olinto 
(1986), several other authors have considered this process as 
a probable mechanism of GRBs (Ma & Xie 1996, Cheng & 



Dai 1996). For example, Ma & Xie (1996) studied the phase 
transition of a neutron star into a "hybrid" star. Cheng & 
Dai (1996) also pointed out that a NS in a low mass X-ray 
binary, after accreting sufficient mass may undergo a sim- 
ilar phase transition and power a weak GRB event. While 
Alcock, Farhi & Olino attempted to explain an electromag- 
netic FB of Qfb ~ 10 45-46 erg, the other authors, crudely 
estimated that the value of Qfb ~ 10 49 erg. Clearly such 
low values of Qfb are highly insufficient to explain the en- 
ergetics of presently observed GRBs. Thus while the basic 
concept involved in this paper is not new (NS — > SS), this 
general relativistic model as such is new because it may ex- 
plain four order higher values of Qfb- 

Since this paper has similarity in idea with that of 
Cheng & Dai (1996), we shall briefly recall the main result of 
the latter paper. Cheng & Dai adopted a different approach 
and estimated that the new born strange star has a thermal 
energy _E t h ~ 5 x 10 52 erg. The star cools by emitting this 
thermal energy in the form of v, v pairs. While the neutri- 
nos traverse out from the core, they impart electromagnetic 
energy inside the star by processes like n + v E — > p + e~ 
and p + v e — > n + e + . Cheng and Dai estimated the optical 
thickness of the above processes to be r ~ 4.5 x 10 -2 and 
concluded that the electomagnetic energy deposited within 
the body of the star to be E2 ~ r x E t h ~ 2 x 10 52 erg , which 
is obviously incorrect for the given value of -Eth and r. Thus, 
as one can easily verify, Cheng & Dai here overestimated E2 
by a factor of 10 and its value works out to be E2 ~ 2 x 10 51 
erg. Moreover, Cheng & Dai overlooked here fact that pairs 
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generated within the body of the star can not come out be- 
cause the star is exceedingly optically thick (electromagnetic 
optical thickness is ~ 2.5 g cm -2 . Thus this electromagnetic 
energy gets absorbed within the star and reemitted primar- 
ily as neutrinos. So the total neutrino flux coming out of the 
star is ~ E t h ~ 5 x 10 52 erg. Now Cheng & Dai esitimated 
that total energy coming out from the star in the form of 
pairs, due to the v + v — > e + + e~ process, is E\ ~ 10 49 erg. 
Thus, the total external pair flux flux energy, in the model 
of Cheng & Dai should be Eo ~ Ej ~ 10 49 erg and not 
Eq = E\ + E-2 ~ 5 x 10 52 erg as erroneously concluded by 
them. 

The general relativistic treatment presented below 
would show that, actually, the pair fireball (FB) enegy can 
be much more. 



2 THE MODEL 

We consider an initial scenario where a NS, either isolated or 
accreting mass in a low mass X-ray binary, is of low compact- 
ness. This is possible for a class of low density stiff NS equa- 
tion of states. For instance, the radius of a 1.4Mq NS with 
a (low) fiducial density parameter po could be (Kalogera & 
Baym 1996). 

-0.35 

PO 



R = 21.2km 



10 14 gcm- 



(1) 



Here po is not to be confused with the central density p c . 
Such low density NSs have a binding energy (BE) insignifi- 
cant compared to their baryonic mass E* B ~ 0.1M Q c 2 . Here 
Mq is solar mass and c is the speed of light. Following Olinto 
(1987) we consider the possibility that the NS is hit by a pri- 
mordial strange matter clump or seed in the form a cosmic 
ray. Once the NS gravitationally captures such a "seed" , the 
neutrons surrounding the seed will be absorbed by it and 
be deconfined to be strange matter (Alcock & Olinto 1988, 
Cheng, Dai & Lu 1998). Thus the strange matter seed can 
become bigger and bigger until the whole star is converted 
into a strange star. And this process may be complete ex- 
tremely fast on a timescale of 10 -7 s (Alcock & Olinto 1988, 
Farhi & Jaffe 1984) This process may occur for an accreting 
NS once it starts collapsing after its mass crosses a certain 
maximum value and the central density increases sufficiently. 
This might also happen for the rapidly spinning supramas- 
sive NS scenario (Vietri & Stella 1998) where the NS starts 
collapsing after loss of sufficient angular momentum. 

The maximum values of masses and radii of strange 
stars may be represented by (Witten 1984, Cheng, Dai & 
Lu 1998) 

M = 2.0M 9 (^) 1/2 (2) 
and 

R = 11-1 (§) V2 ^ (3) 

where the "bag constant" is B and Bo = 57 MeV/fm 3 There- 
fore the surface gravitational redshift of a static spherical 
strange star is given by (Shapiro & Teukolsky 1983) : 



' _ 2GM(R) \ 1/2 
i?c 2 



1 « 0.5 



(4) 



Here r is the invariant circumference radius, c is the speed of 
light, and M is the gravitational mass enclosed within r — r 



M(r) = / pdV = / dM 
Jo Jo 



(5) 



where p is the total mass-energy density, dV = 4-7rr 2 dr 
is coordinate volume element, and the symbol dM is self- 
explanatory. 

The self-gravitational energy of a static relativistic star 
is given by (Shapiro & Teukolsky 1983) 
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Then recalling the definition of z from Eq. (4), we may write 
z{r)dM « azM w -zM (7) 



E. 



where a ~ 1 is a model dependent parameter and studying 
the models of relativistic polytropes, we have numerically 
verified that for z < 0.6, indeed a > 1. 

The binding energy of a cold star is given by Eb ~ 
(1/2) | Eg |. If the NS undergoes a phase transition to be- 
come a more compact SS, the total energy to be liberated 
would be Eb(SS) — Eb(NS). But here we are consider- 
ing only those situations where the original NS is not very 
compact (Kalogera & Baym 1996) with a canonical value 
of Zi < 0.15, and thus, for the sake of analytical simplic- 
ity we neglect the initial B.E. of the NS in comparison to 
the much larger B.E. of the SS, and to compensate for the 
initial B.E. we set here a = 1.0 although, actually it is 
marginally higher. Since occurrence of powerful GRBs hav- 
ing a frequency ~ 10 ~ 6 — 10~ 7 /galaxy /yr (under conditions 
of isotropy) is an extremely rare event compared to the 
same for supernova events, it may be justified to consider 
only favourable initial condtions. Further during the NS- 
SS conversion, an additional energy of ~ 30 MeV/nucleon 
is liberated due to quark deconfinement (Alcock & Olinto 
1988, Fari & Jaffe 1984). For the sake of analytical simplic- 
ity, we club all such liberated energy into a single expression 
Q v ~ Eb ~ *M-. In the case of NS formation, the nascent 
hot NS radiates most of its BE in the form of v — v on a time 
scale of ~ 10 s (Shapiro & Teukolsky 1983) and we take this 
as a fiducial scale for the present problem. 

So, given this value of z = 0.5 the maximum value of 
Q v ta 1.0 x 10 54 M 2 erg where of M = M 2 2M e . The value 
of Q v measured near the compact object will be higher by a 
factor (l + z): Q' v — ^(l + z)M/2. And the locally measured 
duration of the burst would be t! v = (1 + z)~ x t v . Therefore, 
the mean (local) v — v luminosity will be 



L' 



<X = z{\ + z) 2 M ^ 

2 x 10 53 z(l + z) 2 M 2 t-Q 1 erg/s 



(8) 



where t u = tiolOs. As to the cooling of strange stars, there 
is an important difference with respect to a NS: while a NS 
cools via emission of all 3 flavours of u, u of approximately 
equal proportion, a SS cools primarily via the emission of 
electronic v e ,v e : 



u+e 



d+u e , u+e 



s+u e , d — > u+e +v e 



u+e +^ e (9) 
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Thus there will primarily be only two specieses of neutrino 
each contributing to a luminosity L' t — (1/2)L' V . In contrast 
normal baryonic matter (like a NS) cools via emission of 
six specises of neutrinos. Consequently, emissivity of elec- 
tromagnetic pairs is much higher for a SS cooling in com- 
parison to a NS cooling, and we feel that, this point has not 
been noted earlier. By assuming the radius of the neutri- 
nosphere to be R v « R, the value of effective local neutrino 
temperature T' is obtained from the condition 



electromagnetic FB at oo is 



L' = -4-kR 2 <jT" 



(10) 



where a is the Stefan-Boltzman constant. By equating Eqs. 
(8) and (10), we find 



T' = 



' 4z(l + zfMC 



7TvaR 2 t„ 



1/4 
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where R = ReW 6 . For a Fermi-Dirac distribution, under the 
crude assumption of zero i^-chemical potential, the mean 
(local) energy of the neutrinos is E' v w 3.15T' ~ 57 MeV 
(for z = .5). The various neutrinos will collide with their 
respective antiparticles to produce electromagnetic pairs by 
the v-\-u — > e + +e~ process. The rate of energy generation by 
pair production per unit volume per unit time, at a distance 
r from the center of the star, is given by (Goodman, Dar & 
Nussinov 1987) 



<7±M = ^2 



K vl G 2 F E' v L l2 (r) 
12-k 2 cRi 



ip(r) 



(12) 



Here, I4( r ) ~ r 2 is the z^-flux density of a given species 
above the //-sphere, G% = 5.29 x 10~ 44 cm 2 MeV~ 2 is the 
universal Fermi weak coupling constant squared, K V i = 2.34 
for electron neutrinos and has a value of 0.503 for muon and 
tau neutrinos. Here the geometrical factor ip(r) is 

(iV0 2 ] 1/2 (13) 



<fi(r) = (1 -xf(x l +4x + 5); 
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But since in the present case, most of the neutrinos are of 
the electronic nature, there is substantial enhancement of 
the efficiency: 



. , \ K ve G 2 F E v L 2 {r) . . 
9±W = JoZ^Tul ¥>(r) 



(14) 



48n 2 cRi 

Now, a simple numerical integration yields the local value 
of pair luminosity produced above the neutrinosphere : 



L' + 



■A 2 7 K v £ Gf 2 E' V L V 52 

q±VKT dr ~ — — ~ 1.2 x 10 



277raR„ 

:- - (1 + z) AJ > Ml 25 ^o 2 ' 25 Rferg/s (15) 



R 

2.25 



This estimate is obtained by assuming rectilinear propaga- 
tion of neutrinos near the SS. Actually, in the strong grav- 
itational field near the SS surface the neutrino orbits will 
be curved with significant higher effective interaction cross- 
section. Since, most of the interactions take place near the 
//-sphere, for a modest range of z, we may tentatively try 
to incorporate this nonlinear effect by inserting a (1 + z) 2 
factor in the above expression. On the other hand, the value 
of this electromagnetic luminosity measured by a distant 
observer will be smaller by a factor of (1 + z)~ 2 , so that 
eventually, L± — L'± of Eq.(15). And the total energy of the 



Qfb 
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For z = 0.5, M.2 ~ 1, i?6 ~ 1.1, and t = 10s, we obtain a 
highest value of Qfb ~ 1.5 x 10 53 erg. The efficiency for 
conversion of Q v into Qfb in this case is e± = Qfb/Qv ~ 
15%. It is known that NSs have an upper limit of mass 
« 3.OM0 and with the inclusion of rotation the maximum 
mass could be as high as 3.7M0 (Shapiro & Teukolsky 1983). 
Similarly for a bag constant B — Bo — 57 MeV/fm 3 , the 
maximum mass of a rotating SS is 2.4 Mq (Colpi & Miller 
1992). At this stage the actual value of Bo is uncertain and 
thus a value of M > 3Mq can not be ruled out for a SS. Then 
Eq. (16) may yield a maximum value of Qfb ~ 3.7 x 10 53 
erg, as may be required for GRB 971214. 



3 DISCUSSION 

In some ways, our model in the General Relativistic exten- 
sion of the previous work by Cheng & Dai (1996), and it 
is free from the "baryon pollution problem" in same way. 
The mass of the baryonic crust of a canonical SS is only 
AM = 2x 10" 5 M Q (Cheng, Dai & Lu 1998, Alcock & Olinto 
1984). For the massive SS, even if the value of AM is consid- 
erably higher, we would always have the degree of baryonic 
pollution rj = Qfb /AM > 10 2 , and the genesis of a GRB 
may be understood (Meszaros & Rees 1992). 

For quick comprehension by the readers, let us explain 
here, why our relativistic model is able to generate a FB 
energy as high as four orders of magnitude than previous 
non-relativistic estimates. In a non-relativistic model, one 
usually considers a NS with mass of 1 Mq, so that M2 = 0.5. 
The radius of such a NS is taken as 10 Km so that Ra — 1. 
For such a NS z w 0.15. We call this Case I. On the other 
hand, for the present relativistic model (Case II) involving 
a heavy NS with mass 2 Mq, we have M2 = 1. The radius 
of such a NS is 11 Km which is practically the same as in 
Case I, so that for qualitative understanding, we neglect any 
variation with respect to R in the two cases. As mentioned 
before, in this latter case 2 w 0.5. 

Then for a a burst of duration of 10 s, the essential 
relativistic fators in Eq. (16) are 



Qfb ~ 2 2 - 25 (l + 2) 4 ' 5 Mf' 



(17) 



For case I, the numerical value of this factor is Qfb ~ 
0.005. And for Case II, the same numerical factor is Qfb ~ 
1.3. Thus, ignoring higher pair emissivity for the moment, 
the ratio of relativistic factors in the two cases is 260. The 
previous authors ignored the fact that for SS cooling, pair 
emissivity is much higher compared to the case of NS cool- 
ing. For the latter case (Case I), for a total neutrino luminos- 
ity of L' v , the luminosity in each species is L[ = (1/6)L^. In 
contrast for the former case (Case II), we have L\ = (1/2)L^. 
The pair emissivity, as per Eq. (12) is oc L' 2 . The same is also 
<x El txT' <x Lf- 25 . Then finally, pair emissivity oc L' 2 25 . 
Since, for Case II, L[ is higher of 3, one would obtain a pair 
emissivity higher by a factor of w 12 here if the value of K V i 
were same for all the flavours of neutrinos. But, as men- 
tioned before, while K V i — 0.503 for /1 and r neutrinos, its 
value is 2.34 for electronic neutrinos. This fact enhances the 
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pair emissivity approximately by a factor of 4 for Case II. 
Considering all the three effects, FB energy becomes higher 
by a net factor of 260 x 12 x 4 sa 10 4 in Case II! 

It might be possible that upon capturing a sufficiently 
massive strangelet, not only a NS, but an ordinary star like 
the Sun can also get converted into a SS (Alcock & Olinto 
1988). Here the burning process into strange matter arises 
via a series of proton capture reactions apart from the usual 
neutron capture and deconfinement mode (Alcock & Olinto 
1988). If the NS— » SS process is complete within 10~ 7 s, it 
is conceivable that a normal star — > SS process is completed 
well within 1 s. And if a massive star gets converted into a 
maasive SS in this way, we would be able to explain origin 
of GRBs having energy even much higher than 10 54 erg. 
However, the resultant massive SS must be unstable and is 
likely to proceed for further collapse. 

In the case of a supernova, even in an initial spherical 
model there could be symmetry breaking Rayleigh- Taylor 
and Kelvin-Helmoltz instabilities because of the severe den- 
sity gradients in the presence of extreme and varable gravity 
(Burrows 2000). There could be additional asyymetry if the 
SS is magnetized and spins fast, and thus it is plausible that 
the resultant GRB will be beamed. Finally, the recent dis- 
covery of a quark-gluon plasma state in CERN provided as 
an additional motivation behind this work (Abott 2000). 
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